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ABSTRACT 

The magnitude-redshift relation for Type la supernovae is beginning to provide 
strong constraints on the cosmic densities contributed by matter, fim, and a cosmo- 
logical constant, SIa, though the results are highly degenerate in the f2„i-f2A plane. 
Here we estimate the constraints that can be placed on a cosmological constant or 
quintessence-like component by extending supernovae samples to high redshift. Such 
measurements, when combined with constraints from anisotropies in the cosmic mi- 
crowave background, could provide an important consistency check of systematic er- 
rors in the supernovae data. A large campaign of high-z supernovae observations with 
10 metre class telescopes could constraint ^Im to a (Icr) accuracy of 0.06 and f^A to 
0.15. A sample of supernovae at redshift z ~ 3, as might be achievable with a Next 
Generation Space Telescope, could constrain to an accuracy of about 0.02 inde- 
pendently of the value of I^a- The constraints on a more general equation of state, 
wq = vl Pi converge slowly as the redshift of the supernovae data is increased. The 
most promising way of setting accurate constraints on wq is by combining high-z su- 
pernovae and CMB measurements. With feasible measurements it should be possible 
to constrain wq to a precision of about 0.06, if the Universe is assumed to be spatially 
flat. We use the recent supernovae sample of Perlmutter et al. and observations of the 
CMB anisotropies to constraint the equation of state in quintessence-like models via 
a likelihood analysis. The 2cr upper limits are wq < —0.6 if the Universe is assumed 
to be spatially flat, and wq < —0.4 for universes of arbitrary spatial curvature. The 
upper limit derived for a spatially flat Universe is close to the lower limit [wq « —0.7) 
allowed for simple potentials, implying that additional fine tuning may be required to 
construct a viable quintessence model. 



1 INTRODUCTION 

The possible discovery of an accelerating Universe from ob- 
servations of Type la supernovae (Perlmutter et al. 1998, 
Riess et al. 1999) has led to a resurgence of interest in 
the possibility that the Universe is dominated by a cosmo- 
logical constant (for a recent review see Turner 1999). A 
number of authors have shown how observations of distant 
Type la supernovae (SN) can be combined with observations 
of CMB anistropies to constrain the cosmological constant 
and matter density of the Universe (White 1998, Tegmark 
et al. 1998, Lineweaver 1998, Garnavich et al. , 1998, Ef- 
stathiou and Bond 1999, Tegmark 1999, Efstathiou et al. 
1999). For example, Efstathiou et al. (1999, hereafter E99) 
combine the large SN sample of the Supernova Cosmology 
Project (Perlmutter et al. 1998, hereafter P98; we will refer 
to these supernovae as the SCP sample) with a compilation 
of CMB anisotropy measurements and find Q,rn ~ 0.25I^q ]^2 
and = 0.63+0 23 (95% confidence errors) for the cos- 



mic densities contributed by matter and a cosmological con- 
stant respectively. These results are consistent with a num- 
ber of other measurements, including dynamical measure- 
ments of i},n, the large-scale clustering of galaxies and the 
abundances of rich clusters of galaxies (Turner 1999, Bridle 
et al. 1999, Wang et al. 1999). 

The observational evidence for an accelerating Universe 
has stimulated interest in more general models containing a 
component with an arbitrary equation of state, p/p — wq 
with WQ > —1. Examples include a dynamically evolving 
scalar field (see e.g. Ratra and Peebles 1988 and Caldwell, 
Dave and Steinhardt 1998, who have dubbed such a com- 
ponent 'quintessence'; we will refer to this as a 'Q' com- 
ponent hereafter) and a frustrated network of topological 
defect (Spergel and Pen 1997, Bucher and Spergel 1999). In 
particular, Steinhardt , Wang and Zlatev, 1998, have pointed 
out that for a wide class of potentials, the evolution of a Q- 
like scalar field follows 'tracking solutions' , in which the late 
time evolution is almost independent of initial conditions. 
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(a) SCP + 20 SN <z> = l 



(b) 2xSCP + 40 SN <z> = 1.5 




Figure 1. The dashed Unes in each panel show 1, 2 and Sir lilselihood contours in the Q.t^-^m plane for the SCP distant supernova 
sample as analysed by E99. The solid contours are derived from the Fisher matrix (equation 4) for the SCP sample supplemented by 20 
SN with a mean redshift of (z) = 1 (Figure la) and for twice the SCP sample and 40 SN with {z) = 1.5 (Figure lb). The points show 
maximum likelihood values of Ha and Qm for Monte-Carlo realizations of these samples, as described in the text. 



The purpose of this paper is three-fold. Firstly, to illus- 
trate how the constraints on f2m and can be improved 

by extending the redshift range of the supernovae samples. 
At low redshifts, the magnitude-redshift relation is degen- 
erate for models with the same value of the decellaration 
parameter go (= \{^m — 2nA)). This degeneracy can be 
broken by observing supernovae at redshifts ^ 1 (see, for 
example, Goobar and Pcrlmuttcr, 1995). Thus, by extend- 
ing the redshift range of the current supernovae samples it 
should be possible to set tighter limits on fi^ra and Qa in- 
dependently. This is important because there are significant 
worries that the SN data may be affected by grey extinction, 
evolution, or some other systematic effect. The consistency 
of SN constraints on Q,m and Q,k with those derived from the 
CMB anisotropy measurements would provide an important 
consistency check of systematic errors in the SN data and 
the interpretation of the CMB data. Secondly, we estimate 
the accuracy with which a more general Q-likc equation of 
state can be constrained by high redshift SN and CMB data. 
Thirdly, we use the current SN and CMB anisotropy data 
to constrain Q-like models in a spatially flat universe and in 
a universe with arbitrary spatial curvature. 



2 ANALYSIS OF MODELS WITH A 
COSMOLOGICAL CONSTANT 

2.1 Constraints from supernovae at 2: ~ 1 

The predicted peak magnitude-redshift relation is given by 

mF'^'^iz) =M. + 51ogX)z,(2, 0^, Ha), (1) 

where M is related to the peak absolute magnitude by 
M = M - SlogJJo + 25. and Vl = dh + 51oglfo is the 
Hubble constant-free luminosity distance. To compute the 
luminosity distance, we ignore gravitational lensing and use 
the standard expression for a Universe with uniform density 



Table 1: Fisher Matrix Errors, and Ha- 



Supernovae Alone 





SCP 


SCP + 20 SN 


2xSCP + 40SN 


< 2 > 




1.0 


1.5 


1.0 


1.5 




0.53 


0.130 


0.081 


0.092 


0.057 


snA 


0.71 


0.265 


0.218 


0.19 


0.154 


SM 


0.056 


0.053 


0.049 


0.035 


0.035 






Supernovae -|- CMB 






SCP 


SCP + 


20 SN 


2xP98 


f 40SN 


<z> 




1.0 


1.5 


1.0 


1.5 


8Q,jn 


0.073 


0.055 


0.047 


0.039 


0.033 


SQa 


0.080 


0.060 


0.051 


0.042 


0.036 


SM 


0.046 


0.042 


0.039 


0.030 


0.028 



(see e.g. Peebles 1993), 



Jo 



dz' 



[n^(i + z'Y + Q.k{i + z'Y + nA]i/2 



(2) 



where f2fe = 1 — Q.m — f^A and sin^ = sinh if > and 
sink = sin for Q.k < 0. 

We assume that we observe N supernovae, with peak 
magnitude rrii, (corrected for k-tcrm, decline rate- luminosity 
relation, reddening etc), magnitude error ai and redshift Zi, 
from which we want to determine a set of parameters Sk by 
maximising the likelihood function. 



cxp 



prcd\ 2 



20-? 



(3) 



In this section we assume that the parameters Sk are fim, 
Qa and M. (defined in equation 1). An estimate of the co- 
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variance matrix, dj, for these parameters for a given SN 
data set is given by the inverse of the Fisher matrix 



SCP + 20 SN <z>=l 



-I a prcd o pred 

1 dm^ dm); 



^-^ a 

k 



dsi 



dsi 



(4) 



(Kendall and Stewart 1979). The marginalized error on each 
parameter (given by ^/Cu) is listed in Table 1 for several as- 
sumed supernova datascts. The cohnmi labelled SCP gives 
the Fisher matrix errors on Qrn , and M derived for sam- 
ple C (56 supernovae) of P98, i. e. assuming the magnitude 
errors, intrinsic magnitude scatter and rodshift distribution 
of the real sample. The next two columns give the expected 
errors for the SCP sample supplemented by 20 supernovae 
with a peak magnitude error of Am = 0.25 magnitudes and 
a Gaussian redshift distribution of dispersion Az = 0.5 and 
mean redshift (z) — 1 and 1.5. The upper redshift limit 
is close to the maximum for feasible spectroscopic mea- 
surements with 10 metre-class telescopes (see Goobar and 
Perlmutter 1995). As these authors comment, ground based 
spectroscopy at optical wavelengths becomes prohibitively 
expensive for supernovae at higher rcdsliifts because of the 
strong K-correction. The last two columns give the errors for 
a sample twice as large as the SCP sample supplemented by 
40 supernovae with mean redshift of 1.0 and 1.5. Wo adopt 
a background cosmology with SI a = 0.63 and flm = 0.25 as 
indicated by the joint likelihood analysis of the SCP sample 
and CMB anisotropics described in E99. 

Ftom Table 1 we see that the Fisher matrix analysis 
of the SCP sample gives relatively largo errors on ^Ia and 
ilm, in agreement with the likelihood analysis presented by 
P98. However, by adding 20 SN at 2 ~ 1, the errors on 
Oa and in particularly fi^ are reduced significantly. The 
last column shows that an enhanced SCP sample together 
with 40 SN at 2: ~ 1.5 (a formidable, but feasible ob- 
serving programme) can provide a tight constraint on fim. 
The parameters Qa and Qm. are, of course, highly corre- 
lated. This is illustrated in Figure 1 which shows 1, 2 and 
3(T error ellipses in the riA-^m plane after marginalizing 
over S3 — A4 assuming a uniform prior distribution. (The 
components of the new Fisher matrix after marginalization 
are given by Fl^ = Fn - F^s/Fsz, = F22 - FI3/F33, 
FI2 = F12 — F13F23/F33.) The points in the Figure show the 
results of Monte-Carlo calculations, where we have simu- 
lated the observational samples and determined the param- 
eters Si by ma^ximising the likelihood function (equation 2). 
By diagonalizing the matrix F' we can find the orthogonal 
linear combinations fiy = aflm +&0a and n± = 60m — fflf^m 
defining the major and minor axes of the likelihood con- 
tours shown in Figure 1. The distributions in these orthogo- 
nal directions are shown in Figure 2 and compared with the 
distributions determined from the Monte-Carlo simulations. 
The Montc-Carlo distributions are very close to Gaussians 
and show that the Fisher matrix gives an extremely accurate 
description of the errors in the Q.m f^A plane. 

Although the errors in and Qa are significantly re- 
duced by the addition of high redshift supernovae over those 
of the SCP sample, they are still quite large in the parallel 
direction f2|| . This means that it is difficult to set tight limits 
on Qa from SN measurements alone. The constraints on the 
spatial curvature f2fc are even weaJter. For example, for the 
larger sample shown in Figures 1 and 2, the la error on Qk is 




-0.4 -0.2 0.2 0.4 0.5 1 1.5 

n I =o.923n„-o.385n, nii=o.385n„+o.923n, 

± m A II m A 

2xSCP + 40 SN <z> = 1.5 




-0.4 -0.2 0.2 0.4 0.5 1 1.5 

n I =o.967n„-o.254n. n||=o.254n„+o.964n, 

-L m A II m A 

Figure 2. Distributions along the major and minor axes of the 
likelihood contours shown in Figure 1. The histograms show the 
distributions derived from the Monte-Carlo simulations and the 
dotted lines sliow Gaussian distributions with variances deter- 
mined from the Fisher matrix after marginalizing over the pa- 
rameter M. 



SQk =0.19. This can be reduced by extending the range to 
even higher redshifts (see Section 2.2) or by combining the 
SN data with cosmic microwave background anisotropics, as 
has been done by several authors (White 1998, Lineweaver 
1998, Garnavich et al. 1998, Tegmark 1999, E99). 

CMB anisotropy measurements, especially with future 
satellites such as MAP and Planck, are capable of setting 
tight constraints on the locations of the acoustic peaks in the 
CMB power spectrum. Following E99, we define an acous- 
tic peak location parameter 'yoiyimy^A) to be the ratio of 
the peak position in a model with arbitrary cosmology com- 
pared to that in a spatially fiat model with zero cosmological 
constant. (This parameter depends weakly on the matter 
content of the Universe and on the spectral index of the 
fiuctuations, but we ignore these small dependences in what 
follows). CMB measurements are therefore capable of fix- 
ing 7i3, defining a degeneracy direction in the Q.A-^m plane 
given by 
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2XSCP + 100 SN <z>=1.5 + 40 SN <z> = 3 




Figure 3. Fisher matrix constraints for a sample of SN extending 
to redshifts 2 > 3 (see text) illustrating that by extending the 
redshift range one can determine Clm, independently of Q,/^. 



. „ (97£/9an) 
AS2a = -77; /ao N ^^'m, 



(5) 



(see Efstathiou and Bond 1998). The results in the lower 
panel of Table 1 show the Fisher matrix analysis of the SN 
samples including the constraint imposed by equation (5). 
As is well known, the combination of SN and CMB measure- 
ments can break the degeneracy between Q.a and Q,m and 
it should be possible to determine these parameters with an 
error of less than 0.04 with an enlarged supernova sample 
assuming, of course, that systematic errors are unimportant. 

Although the errors on Q,a from SN measurements alone 
converge relatively slowly as the redshift range is increased, 
consistency of the cosmological parameter estimates pro- 
vides a strong test of systematic errors in the SN data. If 
we believe that systematic errors are unimportant, and that 
our interpretation of the CMB anisotropies (in terms of adi- 
abatic CDM-liko models) arc correct, then current data al- 
ready constrain Q.^ and JIa to high precision (see Fig 5 of 
E99). Consistency requires that the likelihood contours for a 
high redshift supernova sample converge to the same answer. 



2.2 Constraining 0,„ with NGST 

Observations of very distant supernovae a± z ^'i may be 
possible with a Next Generation Space Telescope {e.g. 
Miralda-Escude and Rees 1997, Madau 1998, Livio 1999). 
We will not analyse the feasibility of such observations here. 
Rather, we note from Figures 1 and 2 that the major axis of 
the error ellipses in the JIa- Q.m tilt and become more verti- 
cal as the redshift range of the SN sample is increased. This 
is because the magnitude redshift relations for models with 
very different values values of Qa and the same Om con- 
verge at higher redshifts. The convergence redshift depends 
on Q,m and lies between z « 2-4 for Q,m in the range 0.2-1 
(see Figure 1 of Melnick, Terlevich and Terlevich, 1999). 

This is illustrated by Figure 3, which shows the 1, 2 
and 3(7 likelihood contours determined from the Fisher mar 



trix for a sample consisting of twice the SCP sample, 100 
SN with {z) = 1.5, Az = 0.5, and 40 SN with (2) = 3, 
= 1. As expected, these contours are almost vertical in 
the Q,m-^h plane. A sample of supernovae (or some other 
distance indicator such as HII galEixies, Melnick et al. 1999) 
at redshifts a: ~ 3 can therefore produce a tight constraint 
on Q,Tn independently of the value of f^A. 



3 CONSTRAINTS ON AN ARBITRARY 
EQUATION OF STATE 

In this Section, we analyse the constraints that SN can place 
on an arbitrary equation of state. We first consider a con- 
stant equation of state. Models of this typo (sec Bucher 
and Spergel) include a frustrated network of cosmic strings 
{p/p = —1/3) and a frustrated network of domain walls 
(p/p = —2/3). A constant equation of state is also a good 
approximation to a Q component obeying tracker solutions. 
Tracker solutions arc discussed in Section 3.2. Constraints 
on generalised forms of dark matter with anisotropic stress 
are discussed by Hu et al. (1999) and will not be considered 
here. 



3.1 Constant equation of state 

If we include a Q-like component with equation of state 
p/p = WQ, the expression for the term x in the luminos- 
ity distance (equation 2) is modified to 



dz' 



[n„(l + z')'' + f^fc(l + ^')' + f^Q(l + ^')'''+"''3)]i/2- 

The addition of the parameter wq means that it is not 
possible to constrain all of the parameters fim, S^q, wq to 
high accuracy from the supernova data alone (see Section 
4.2). Thus, Garnavich et al. (1998) analyse the High-z Su- 
pernovae Search (HZS) sample (Riess et al. 1999) assuming 
a spatially flat universe and find that wq < —0.55 at 95% 
confidence. A similar analysis of the SCP sample by Perl- 
mutter, Turner and White (1999) yields wq ^ — 0.5. 

Table 2 lists the results of a Fisher matrix analysis for a 
Q-like component with a constant wq . Here we have applied 
the constraints wq > — 1 and flQ > 0. The upper table 
gives results for the supernovae magnitude-redshift relation 
alone assuming a spatially flat Universe with Clm = 0.25 and 
Wq = — 1. The constraints on wq from a sample such as the 
SCP data are quite poor and improve relatively slowly as 
the sample is extended to higher redshift because of a strong 
degeneracy between wq and flm in the magnitude-redshift 
relation. This is illustrated in Figure 4, which shows the 
analogue of Figure 2 for Q-like models. As the supernovae 
sample is extcuded to higher redshift, the likelihood contours 
narrow but wq and ^Im remain strongly degenerate. 

The situation is dramatically improved by the addition 
of constraints from CMB anisotropies. The addition of a 
Q-like component affects the location of the Doppler peaks 
(see Caldwell et al. 1998, White 1998) and, in analogy with 
equation (5) , an accurate determination of the CMB power 
spectrum imposes the constraint 
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(a) SCP + 20 SN <z>=l (b) 2xSCP + 40 SN <z> = 1.5 





Figure 4. As Figure 2, but for an arbitrary constant equation of state in a spatially flat Universe. The dashed lines in each panel show 
1, 2 and 3a likelihood contours in the WQ-Qm plane for the SCP distant supernova sample as analysed in Section 4 (assuming a constant 
equation of state) . The solid contours are derived from the Fisher matrix for enhanced samples of high redshift supernovae and the points 
show maximum likelihood derived from Monte-Carlo realizations of these samples. 



Table 2: Fisher Matrix Errors, Om, and wq. 



Supernovae Alone, Qfc = 





SCP 


SCP + 20 SN 


2xSCP + 40SN 


< z> 




1.0 


1.5 


1.0 


1.5 


Sflm 


0.14 


0.12 


0.097 


0.097 


0.073 


5WQ 


0.36 


0.35 


0.32 


0.28 


0.24 


SM 


0.051 


0.051 


0.051 


0.037 


0.036 




Supernovae 


f CMB, Qk 


= 






SCP 


SCP ^ 


- 20 SN 


2xSCP 


+ 40SN 


<z> 




1.0 


1.5 


1.0 


1.5 




0.027 


0.022 


0.0210 


0.016 


0.015 


SWQ 


0.10 


0.085 


0.081 


0.061 


0.057 


SM 


0.048 


0.046 


0.045 


0.032 


0.032 




Supernovae - 


h CMB, Qf 


= 5^0 






SCP 


SCP + 20 SN 


2xSCP 


+ 40SN 


<z> 




1.0 


1.5 


1.0 


1.5 




0.14 


0.12 


0.095 


0.094 


0.069 


SVLq 


0.10 


0.083 


0.066 


0.067 


0.048 


SwQ 


0.31 


0.31 


0.27 


0.24 


0.20 


SM 


0.051 


0.051 


0.051 


0.036 


0.036 



(7) 



The second panel of Table 2 shows the constraints derived 
on an arbitrary equation of state by combining supernovae 
data with the CMB constraint of equation (7). For spatially 
flat models, the combination of SN and CMB anisotropics 
constrains wq to an accuracy of better than 0.1, sufficient 
to set tight constraints on the physical parameters of Q- 
like models (for example, whether one requires contrived 
potentials, see Section 4). However, the constraints on wq 
improve relatively slowly as the SN sample is extended to 
higher redshift. Similar conclusions apply if the assumption 
of a spatially fiat universe is relaxed (see the lower panel 



of Table 2). In that case, the parameters fl-m and fig can 
be determined to high precision, but the constraints on wq 
improve slowly as the SN sample is increased. This imphes 
that it is worth analysing the constraints on Q-like models 
with arbitrary spatial curvature using current SN and CMB 
data (see Section 4.2). 



3.2 Time varying equation of state: tracker 
solutions 

In the previous section we have investigated the simplified 
case of a constant wq. If, in fact, the Q-like component 
arises from a slowly rolling scalar field evolving in a potential 
V{Q), the equation of state of the Q component will vary as 
a function of time. The equations of motion of the Q field 
can be written in the following compact form (Steinhardt, 
Wang and Zlatev, 1998) 



V"V 



1 + 



Wb — WQ 

2(1 + wq) 



1 + wb — 2wQ X 
2(1 + wq) 6 + x 

_ (I + wq) 



(8) 



(1 + wq) (6 + a;)2' {1 - wq) 

where primes denote derivatives with respect to Q, a; = 
dlnx/dlna, x — (fhix/dhi^a and a is the scale factor of the 
cosmological model. For a wide cIeiss of potentials, and al- 
most independently of the initial conditions, the evolution 
of Q locks on to a tracking solution in which Q and wq vary 
slowly (see Zlatev, Wang and Steinhardt 1998, Steinhardt et 
al. 1998). Examples of the evolution of wq and Qq at late 
times are shown in Figure 5 for three forms of the poten- 
tial V{Q). In each case, the evolution of wq at ^ 4 is well 
approximated by 



wq = wq(oo) + aln{a/ao) 



(9) 



where a is a small number determined from the value of x 
at the present time. 
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Figure 5. The evolution of the equation of state wq and its 
contribution to the cosmic density parameter fig as a function 
of redshift derived from solutions to the tracker equation (8) for 
three potentials: V{Q) = M''(exp(l/Q) — 1) (figures 5a and 5b); 
V(Q) = W^iM/Qf (figures 5c and 5d); V{Q) = M'^(M/Qf 
(figures 5e and 5f). The curves in each figure are computed by 
varying the parameter M, with more negative values of wq cor- 
responding to higher values of Q.q. 

Figure 7 shows the relations between Q.q, wq and a at 
the present time derived from the solutions to equation 4 
for the three potentials considered in Figure 5. The mini- 
mum value of a is about —0.14, reflecting the fact that Q is 
evolving relatively slowly even at late times. 

With the approximation of equation (9), the energy 
density of the Q component evolves according to 

— T-T = {—) exp --a In a/ao . 10 

Note also that with the approximation of equation (9), the 
tracker equation (8) becomes an algebraic equation relating 
Jyi^'A to WQ, Q,Q and a {wb = wqQq in the matter domi- 
nated era). 

A small value of a ~ —0.1 to —0.2 cannot be deter- 
mined accurately from SN and CMB observations because 
it is highly degenerate with wq and flm- As we will show 



Figure 6. The left hand panels show the tracker solution relations 
between Qq and wq at the present day for the three potentials 
used in Figure 5. The right hand panels show the derivative a = 
dwQ /dlna for the tracker solutions as a function of wq . 

in the next Section, the introduction of the parameter a 
provides a convenient way of testing the sensitivity of con- 
straints on Q-like models to the time evolution of wq. 

We note that Huterer and Turner (1998) have recently 
proposed a prescription for reconstructing the potential of a 
Q-like component directly from the magnitude-redshift re- 
lation of Type la supernovac. This approach may produce 
interesting constraints if the field Q is rapidly evolving at 
late times. For tracker solutions, however, the equation of 
state changes so slowly that it would be difficult to distin- 
guish the true potential from a perfectly flat one. 



4 LIMITS ON THE EQUATION OF STATE 
FROM TYPE lA SUPERNOVAE AND THE 
COSMIC MICROWAVE BACKGROUND 

4.1 Spatially flat models 

In this Section, we use current SN and CMB data to con- 
strain the equation of state of the Universe. The analysis 
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0.2 ' 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



m m 

Figure 7. Constraints of wq and f2m for spatially fiat universes. Figure 7a sliows results for the SCP supernova sample following a 
similar analysis to that presented by E99. Figure 7b shows results for the supernova sample combined with the constraints from the 
CMB anisotropy measurements as described in E99. The contours show 1, 2 and 3a" likelihood contours. The solid contours are derived 
for a = 0, dotted contours are for a = —0.1 and dashed contours for a = —0.2. 



closely follows that presented in E99. We use the sample 
of 56 Type la SN of fit C of P98 and adopt the likelihood 
analysis described by E99 (including a parametric fit to the 
luminosity-decline rate correlation), modifying the expres- 
sion for luminosity distance to incorporate the parameters 
of the Q-like model. The CMB data that we use are plot- 
ted in Figure 1 of E99. We perform a likelihood analysis for 
these data assuming scalar adiabatic perturbations, varying 
the amplitude of the fluctuation spectrum, the scalar spec- 
tral index, the physical densities of the CDM and baryons 
= fic/i^, ujh = Q.bh? ^ and the Doppler peak location pa- 
rameter 7_D. Modiflcations to the CMB power spectrum aris- 
ing from spatial fluctuations in the Q component are ignored 
as these are negligible in the slowly evolving Q models con- 
sidered here (see Caldwell et al. 1998, Huey et al. 1998). We 
integrate over the CMB likelihood assuming uniform prior 
distributions of the parameters to compute a marginalized 
likelihood for 7d as described in E99. The likelihood func- 
tions for the parameters wq, Q.m and Q.q presented below 
are constructed from the expression for the angular diame- 
ter distance to the last scattering surface and the probability 
distribution of 7d- 

Figure 7 shows the constraints on wq and Q.m for spa- 
tially flat universes. The different line types show the con- 
straints for three different values of the parameter a char- 
acterising the evolution of wq, a = G (solid lines), a = —0.1 
(dotted lines) and a = —0.2 (dashed lines). As described in 
the previous section, these values span the range found for 
tracker solutions for a variety of potentials. These rates of 
evolution are so low that they have very little effect on the 
likelihood contours. The constraints plotted in Figure 7 are 
in very good agreement with those derived by Garnavich et 
al. (1998) from an analysis of the HZS sample, and with the 
analysis of the SCP sample (Perlmutter, Turner and White 

* h is the Hubble constant in units of 100 kms~^Mpc~^. 



1999) and of the combined HZS and SCP samples (Wang et 
al. 1999). The fact that the constraints are weakly depen- 
dent on the size of the SN sample is a consequence of the 
strong degeneracy between wq and Q.m discussed in Section 
3.2. 

Figure 7b shows the results of combining the SN likeli- 
hoods with those determined from the CMB. The likelihood 
peaks at wq = —1, fim = 0.29. Qualitatively, these results 
are similar to those of Perlmutter et al. (1999); the favoured 
cosmology has an equation of state wq = — 1 and wq is con- 
trained to be less than —0.6 at the 2a level. However, in de- 
tail, the constraints in Figure 7b are somewhat less stringent 
than those of Perlmutter et al.,, allowing a broader range in 
Q.m (0.15 ^ Q.m ~ 0.5 at the 2a level). This is because Perl- 
mutter et al. include constraints on the power spectrum of 
galaxy clustering based on the data compiled by Peacock 
and Dodds (1994)[|. In our view this is dangerous because 
it requires a specific assumption concerning the distribu- 
tion of galaxies relative to the mass. Qualitatively, for nearly 
scale-invariant adiabatic models, galaxy clustering imposes 
a constraint on the parameter combination F = Q.mh of 
0.2 ^ F 0.3, if galaxies are assumed to trace the mass fluc- 
tuations on large scales (Efstathiou, Bond and White 1992, 
Maddox, Efstathiou and Sutherland 1996). Combined with 
measurements of the Hubble constant (for which Perlmut- 
ter et al. adopt h ~ 0.65 ± 0.05), galaxy clustering leads 
to a constraint of 0.25 ^ Om ^ 0.5, partly breaking the de- 
generacy between wq and Q.m. The combined SN and CMB 
analysis in Figure 7b provides constraints which are nearly 
as tight, but are much less model dependent. 

The constraints of Figure 7b place strong limits on Q- 
like models. For tracking solutions, the constraint wq ^ 



t Perlmutter et al. (1999) do not combine the SN and CMB likeli- 
hoods but analyse the SN data assuming a spatially flat Universe. 
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Figure 8. Analogue of Figure 5, but for quintessence models with arbitrary spatial curvature. Figures 7a and 7b show marginalized 
likelihoods in the WQ-Clm and Q.Q-Qm planes derived from Type la supernova*. Figures 7c and 7d show the combined likelihoods for 
the Type la and CMB anisotropics. As in Figure 5, the solid contours are derived for a = and dashed contours for a = —0.2. 



—0.6 excludes steep potentials {e.g. V{Q) oc Q^'' with 
(3^2) and the data clearly favour a standard cosmologi- 
cal term [wq = — 1). These limits on wq are very close to 
the lower limit (wq ^ — 0.7) allowed for 'physically well 
motivated' tracker solutions (Steinhardt, Wang and Zlatev, 
1998, i.e. smooth potentials with simple functional forms). 
With a slight improvement of the observations one may be 
forced to fine-tune the shape of the potential to construct a 
viable quintessence model. 

The constraints of Figure 7b are somewhat stronger 
than those of Wang et al. (1999), who perform a 'concor- 
dance analysis' of Q-like models using a number of observa- 
tional constraints including those from Type la supernovae 
and CMB anisotropics. These authors conclude limits of 
— 1 ,$ WQ ^ — 0.4. The difference is caused by the different 
methods of statistical analysis. The concordance analysis of 
Wang et al. leads, by construction, to more conservative lim- 
its than the maximum likelihood analysis and is more robust 
to systematic errors in any particular data set. However, pro- 
vided systematic errors axe negligible in the CMB and SN 
datasets, then the constraints of Figure 7b derived by com- 



bining likelihoods should be realistic. These small differences 
in the upper limits on wq are important because they can 
place significant restrictions on the physics. As stressed in 

the previous paragraph, the upper limit of wq ~ —0.6 places 
strong constraints on tracker models with simple potentials. 

4.2 Models with arbitrary spatial curvature 

Figure 8 shows the results of a likelihood analysis of the 
SN and CMB data, but now allowing arbitrary spatial cur- 
vature. We show two projections of the likelihood distribu- 
tions, marginalizing over ^Iq in Figures 8(a) and 8(c) and 
over Wq in Figures 8(b) and 8(c). The constraints, although 
weaker than those presented in Figure 7, are interesting nev- 
ertheless. The combined SN and CMB likelihoods give a 2cr 
upper limit on wq of wq ^ — 0.4 (Figure 8d) and a maxi- 
mum likelihood solution of flm = 0.12, Q.q = 0.73 irrespec- 
tive of the value of wq. Evidently, the SN and CMB data 
constrain us to a nearly spatially flat Universe dominated 
either by a cosmological constant, or a Q-like component 
with an equation of state wq ^ — 0.4. 



© 0000 RAS, MNRAS 000, 000-000 



hiquanon oj state jrom bupernovae and ulVW Ubservations y 



5 CONCLUSIONS 

Observations of distant Type la supernovae have provided 
important evidence that the Universe may be dominated 
by a cosmological constant (P98, Riess et al. 1999). How- 
ever, the constraints in the fiA-fim plane from current data 
are degenerate along a line defined by Q.a ~ 0.32 + 1.43f2m 
(Figure 1). This degeneracy can be reduced significantly by 
extending the redshift range of the supernovae sample. For 
example, with 20 additional supernovae at redshift z ~ 1.5 
the errdtrs iu Vim tuid SIa could by rydutiyd to 6Q,m ~ 0.08 
and SQ. \ pj 0.22. A sample of supernovae Ai ^ ^ 3 could pru- 
vide an accurate estimate of Qm that is independent of the 
value of f^A. 

The combination of supernovae and CMB anisotropy 
measurements can break the degeneracy between Qa and 
Qm if the initial fluctuations are assumed to be adiabatic 
and characterised by a smooth fluctuation spectrum. This 
method applied to recent supernovae and CMB data sug- 
gests a nearly spatially flat universe dominated by a cosmo- 
logical term with SIa ~ 0.65 (Lineweaver 1998, Garnavich 
et al. 1998, Tegmark 1999, E99). The only plausible way of 
avoiding this conclusion is to appeal to some systematic ef- 
fect in the supernovae data, for example, grey dust or an 
evolutionary effect in the supernovae data such as a metal- 
hcity dependence (see e.g. P98 for a discussion). The degen- 
eracy breaking afforded by extending the supernovae data 
to higher redshift would provide an important consistency 
check of such systematic effects and also on the interpreta- 
tion of the CMB anisotropy data 

The constraints on quintessence-like models with an 
equation of state wq = p/p improve relatively slowly as the 
supernovae data are extended to higher redshift. The most 
promising way of constraining wq seems to be to combine 
supernovae and CMB measurements. We have carried out 
a joint likelihood analysis of CMB anisotropy observations 
and the SCP supernovae data. For a spatially flat Universe 
we derive a 2a upper limit of wq — —0.6. For universes of 
arbitrary spatial curvature, the 2cr upper limit is wq — —0.4. 
The combined SN and CMB likeHhood peaks at ilm = 0.12 
and Qq = 0.73 irrespective of the value of wq, suggesting 
that the Universe is almost spatially flat. The 2a upper limit 
of WQ — —0.6 for spatially flat Universes is close to the min- 
imum value of WQ « — 0.7 allowed for simple quintessence- 
models. This suggests that some fine tuning of the potential 
may be required to construct a viable quintessence model. 
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